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Abstract We calculate the magnetic susceptibil- 
ity in the pseudogap phase of cuprates using the two- 
dimensional t-J model and the Mori projection op- 
erator technique. In this phase, the Fermi arcs lead 
to a quasi-elastic incommensurate magnetic response 
for low temperatures. The theory accounts for the re- 
orientation of the susceptibility maxima from the ax- 
ial to the diagonal direction occurring at small hole 
concentrations in I^-zSr^CuO^ A small cusp of 
the hole dispersion near the Fermi level, which is 
connected with the spin-polaron band, affords the 
growth of the maxima with decreasing frequency. As 
in the superconducting phase, the susceptibility max- 
ima have an hourglass dispersion. The assumption of 
the phase separation into regions of superconducting 
and pseudogap phases in underdoped cuprates ex- 
plains the lack of the superconducting gap in their 
susceptibility. The observed rapid change of widths 
of the susceptibility maxima with frequency is related 
to the contribution of the pseudogap phase. This con- 
tribution may also explain the strengthening of the 
quasi-elastic response by modest magnetic fields and 
by small impurity concentrations. 

1 Introduction 

The magnetic susceptibility of underdoped cuprates 
differs from that of more heavily doped crystals by 
the lack of the so-called superconducting gap - the 
suppression of the imaginary part of the susceptibility 
X" at low frequencies u) in the superconducting state 



[HOG]. In underdoped samples, for temperatures 
below the superconducting transition temperature T c 
the susceptibility is finite up to the lowest measured 
frequencies, and in some experiments it demonstrates 
even the tendency for a growth with decreasing u> [2] . 
This behavior is reflected in the quasi-elastic mag- 
netic neutron scattering observed in lanthanum and 
yttrium cuprates. 

In the considered range of hole concentrations x 
and temperatures T the crystals are either in the su- 
perconducting or in the pseudogap phases. The su- 
perconducting gap in the susceptibility is well repro- 
duced in different theoretical approaches (see, e.g., 
[H El IS]). To our knowledge the mechanism of the 
quasi-elastic magnetic response was much less stud- 
ied. In Ref. [2], based on the sharp difference in 
the spin correlation lengths below and above some 
small frequency it was supposed that below T c the 
underdoped crystal La2- x Sr x Cu04 (LSCO) is phase- 
separated into superconducting and non-supercon- 
ducting regions. The latter regions are in the pseu- 
dogap phase and it can be inferred that they are re- 
sponsible for the low-frequency magnetic response. 

The mechanism of this response may be the follow- 
ing. In comparison with the superconducting phase 
the pseudogap state is characterized by an increased 
damping of hole states and by Fermi arcs - frag- 
ments of the Fermi surface near the nodal points 
[8] . For frequencies below the resonance frequency 
co r the locations of maxima of the magnetic suscepti- 
bility in the Brillouin zone are determined by peaks 
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of the spin-excitation damping [5]. The main con- 
tribution to these peaks is made by the process of 
spin-excitation transformation into a fermion pair. In 
the superconducting phase the respective states are 
located on equal energy contours of the fermion dis- 
persion, the size of which decreases with decreasing 
•jj. This leads to a susceptibility, which vanishes be- 
low some frequency. In contrast, in the pseudogap 
phase due to the Fermi arcs the length of the equal 
energy contours remains finite up to zero frequency, 
that ensures a finite susceptibility up to cu ~ T. Thus, 
indeed it can be expected that in the phase-separated 
state the low-frequency susceptibility is determined 
by regions of the pseudogap phase. 

This picture permits to reach some other conclu- 
sions about the magnetic response of the pseudogap 
state. In the superconducting phase, for frequencies 
immediately below oj r the number of fermion pairs, 
which contribute to the susceptibility on the edges of 
the Brillouin zone, is larger than that contributing on 
the diagonals [3] ■ As a consequence the maxima 
of x" are located on the zone edges at the momenta 
(ir ± k, 7r), (7t, 7r ± k), as it is indeed observed in mod- 
erately doped cuprates [IQ1III] (here and hereafter we 
set the lattice period a of a Cu02 plane as the unit 
of length). The same situation occurs in the pseudo- 
gap phase, in which for moderate x the susceptibility 
maxima are peaked at the zone edges down to lj = 0. 
However, in LSCO the Fermi arc length is approxi- 
mately proportional to the hole concentration [8] . As 
the arc length decreases, the number of fermion states 
contributing to the susceptibility on the zone edges is 
diminished faster than that contributing on the diag- 
onals. As a result for small x and uj the susceptibility 
maxima shift to the diagonals of the Brillouin zone, 
as it is indeed observed in the pseudogap phase of 
LSCO for x < 0.05 [12] ■ For these small hole concen- 
trations as well as for larger x the dispersion of the 
maxima has an hourglass shape, in agreement with 
experimental observations [12) . 

If the hole dispersion has no peculiarities at the 
Fermi level, the susceptibility \" m the pseudogap 
phase decreases monotonously with frequency, re- 
maining finite up to w ~ T. In contrast, in Ref. [2] a 
growth of the susceptibility of an underdoped LSCO 
was observed as the frequency decreases. In the con- 



sidered theory such a behavior can be obtained if 
there is a cusp in the hole dispersion at the Fermi 
level. Such a cusp was found in calculations carried 
out in the two-dimensional (2D) t-J model at moder- 
ate x [13] ■ In this case the cusp is connected with the 
spin-polaron band. A similar cusp is supposedly ob- 
served in some photoemission spectra of underdoped 
cuprates (see, e.g., Fig. 1 in Ref. [14]). 

In the pseudogap phase, the Fermi arcs produce 
a peculiar shape of the susceptibility \'\ which as 
a function of wave vector has a maximum and an 
extended shoulder. The shoulder is weak for small 
frequencies. However, it grows rapidly with u) and 
determines the width of the spectrum for uj > 5 meV. 
As a consequence the width changes steeply near this 
frequency, in close analogy with experimental obser- 
vations in LSCO [2]. 

The structure of the article is the following. The 
main formulas used in the calculations are given in 
the next section. Obtained results are considered in 
Sec. [3] Section [4] is devoted to the discussion of pos- 
sible applications of the obtained results for the in- 
terpretation of the strengthening of the quasi-elastic 
response by the magnetic field and impurities, and to 
the concluding remarks. 



2 Main formulas 

General formulas for the magnetic susceptibility of 
the 2D t-J model in the superconducting phase were 
derived in Ref. [S] with the use of the Mori projection 
operator formalism [15) . The reformulation for the 
pseudogap phase needs only minor changes. There- 
fore, we adduce merely equations, which were used in 
the calculations. 

In the Mori formalism the susceptibility is de- 
scribed by the equation 

X(kw) = - — ^— j . (1) 

In this equation, 

&k = 4(iFi + JCi)(<yk-l), (2) 
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where k is the 2D wave vector, t and J are the Z( kk ') 
hopping and exchange constants between neighbor- 
ing sites, 7k = [cos(fc x ) + cos(fc y )] /2, 
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are correlators of the hole aj^ and spin-i oper- 
ators on the neighboring sites, a = ±1 is the spin 
projection, the angular brackets denote the statistical 
averaging and N is the number of sites. The param- 
eter a serves for correcting the decoupling procedure 
[161 117) . For small x it is approximately equal to 1.7 
[15] , The parameter S describes the gap in the spin- 
excitation spectrum near the antiferromagnetic wave 
vector Q = (tt, tt) due to the finite temperature [17] 
and/or the hole concentration [TB]. The magnitude 
of this gap is directly connected with the correlation 
length of the short-range antiferromagnetic order and 
it grows with x and T. In some cases the magnitude 
of the gap determines the frequency uj t of the waist in 
the hourglass dispersion of the susceptibility maxima 
i- 

There are four contributions to the polarization op- 
erator n(kw) of spin excitations in Eq. (JTJ. These 
contributions are connected with the processes of the 
spin-excitation decay into an electron-hole pair, the 
decays into a pair assisted by a hole or by another 
spin excitation, and the decay into three spin excita- 
tions. For small frequencies considered in this article 
the first of the mentioned processes dominates. Its 
contribution to the imaginary part of the polariza- 
tion operator reads 



t' is the hopping constant between next-nearest 
neighbor sites, 7 k = cos(fc x ) cos(fc a ), i k — 4^7k + 
4i'7 k , N f (uj,uj') = n F (uj) - n F (uj'), n F (uj) = 

[& ) I T + l) . Equation ((4]) was obtained in the sup- 
position that the hole spectral function is described 
by the expression 
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where 



1 if e Q - 2t < oj < -A, 
S(w) = { or A < u> < eq + 6i, 
in other cases, 
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e k is the hole dispersion and A is the pseudogap 
magnitude. The first term on the right-hand side of 
the spectral function ([7]) corresponds to the coherent 
peak and the second term to the incoherent contin- 
uum of the hole-spin-excitation scattering. Such a 
form of the A (lew) is typical for the t-J model (see, 
e.g., |18j). The continuum spans the frequency range 
approximately equal to the width of the initial band 
8t. For moderate doping the lower edge of the contin- 
uum lies approximately 2t below the bottom of the 
hole band £k, which is located at the momentum Q 
in the used hole picture. The spectral weights of the 
coherent and incoherent parts of the spectral function 
are connected by the relation [9] 



Z + Z' = 



1 



To simplify further calculations we set T — (we 
shall discuss possible consequences of a nonzero tem- 
perature in the next section). For T — and small 
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Figure 1: The hole dispersion in the pseudogap phase, 
Eq. (fTOj) . in the first quadrant of the Brillouin zone 
for I = 1 and A = O.lt. Vertical fragments of the 
dispersion near the Fermi level correspond to gapped 
regions. 

positive frequencies the quantity 

/oo 
Nf(u>' 'j OJ + qj')S(lj )S(oj + uj')duj' 
-oo 

in Eq. © reads 

f ifw<2A, 
[UJ> \ lj-2A if 2A < uj < 2t - eq + A 

The real part of n(kcj) is calculated using the 
Kramers-Kronig relation. 

To imitate the hole dispersion in the pseudogap 
phase we use the band 
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cos(k x )+cos(k y ) — 0.3cos(fca;) cos(fcj,)— 0.2, 

(9) 

which has the Fermi surface similar to that observed 
in Lai.88Sro.i2Cu04 [SJ, and construct e k as 



£k = 



if E k > 0, 



Al if£ k <0. 



(10) 



In Eq. (|10p . Ak has the shape of the d-wave gap 
function, Ak = A[cos(k x ) — cos(fc J/ )]/2, everywhere 
except circular regions around the nodal points q = 
(±/t, ±re), Eq — 0. In these regions Ak = 0. Appar- 
ently the diameter I of these regions is close to the 
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Figure 2: Momentum dependencies of x"(kw) for uo = 
O.OOli along the edge [blue solid lines, k = (it, «)] and 
along the diagonal [red dashed lines, k = (ft, «)] of the 
Brillouin zone. The parameter I = 2 in panel (a) and 
0.6 in panel (b). Other parameters are given in the 
text. 



arc length. The hole dispersion £k obtained in such 
a way is shown in Fig. [T] 



3 Results and discussion 

Typical momentum dependencies of the imaginary 
part of the susceptibility obtained with the above 
formulas for a small frequency are given in Fig. [5] 
For this figure the following parameters were used: 
J = 0.2t, 6 = 0.25, A = 0.05i and Z/Z' = 0.1. 
The parameter S was chosen to obtain the frequency 
of the hourglass waist ui r w 40 meV, which is close 
to the value observed in moderately doped cupratcs 
[12] , This frequency is near the magnitude of the 
gap in the spin-excitation spectrum at Q, mentioned 
in the previous section, if u r falls into the region of 
small spin-excitation damping below the edge of the 
electron- hole continuum at Q [5] ■ Here and hereafter 
changing to energy units we proceed from the esti- 
mate t = 500 meV HUD]. Thus : J = 100 meV and 
A = 25 meV, which are also close to the values in 
moderately doped cuprates [5J[T2]. Ratios Z/Z' ~ 0.1 
were obtained in the self-consistent calculations in the 
2D t-J model for hole concentrations 0.07 5, x < 0.15 
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[2"T] . In LSCO the Fermi arc length decreases nearly 
linearly with x 8 . Using the experimental data of 
Ref. [5] we can relate the value I = 2 to x w 0.12 and 
I = 0.6 to x ~ 0.04. For the latter hole concentration 
the parameter 8 has to be decreased, since the value 
of the resonance frequency decreases with x , 
while the pseudogap magnitude A has to be increased 
[5J[2I]. However, for small frequencies these changes 
lead only to minor modifications in the dependencies 
in Fig. [2J 

In the calculations, we have approximated the 5- 
function in the first term of Eq. (|4]) by the step func- 
tion 

[ m other cases. 

It is one of the representations of the <5-function, 
which describes the convolution of the coherent peaks 
in A(kuj) better than the usually used Lorentzian rep- 
resentation, if the peaks have no extended tails. As 
in the superconducting phase, if the Lorentzian rep- 
resentation is used, the incommensurate magnetic re- 
sponse is obtained only for a pronounced nesting of 
the Fermi surface. The LI-shaped representation al- 
lows one to get such a response without the nesting. 
This result is in conformity with experimental obser- 
vations of the incommensurate response in crystals, 
in which Fermi surfaces derived from photoemission 
have no obvious nesting (see Ref. [5] and references 
therein). As follows from Eq. ([5]), we set the intensity 
of the incoherent continuum equal to zero in the range 
— A < io < A. An increase of this intensity leads 
to a growth of x"(Q), which finally suppresses the 
incommensurate peaks when the intensity becomes 
comparable to that outside the range. In the calcula- 
tions, we set rj = ui/3 for uj < 0.005i and rj = 0.002£ 
for larger frequencies. An increase of the width r\ as 
well as the Fermi arc length results in the growth of 

x"(Q). 

As seen in Fig. [5J the low-frequency magnetic re- 
sponse is incommensurate - the maxima on the sym- 
metry lines of the Brillouin zone are displaced from 
Q. For the considered small frequency this incom- 
mensurability is connected with maxima in the spin- 
excitation damping Q descended from the first term 
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Figure 3: The contour plot of the dispersion £k for 
/ = 2 and A — 0.05£. Contours correspond to the 
energies ±0.03t. The squares shown by thin solid 
and dotted lines are the first and magnetic Brillouin 
zones. The blue dashed and red dash-dotted lines 
connect electron-hole pairs contributing to the sus- 
ceptibility on the edge and on the diagonal of the 
Brillouin zone, respectively. 

on its right-hand side. Indeed, from Eq. ([TJ we get 

"fk 1 - a;/iklmn(kii;) 

A ' ~ ' [w 2 - wRen(kw) -iolf + [wlmn(ka;)] 2 ' 

(11) 

The susceptibility above the resonance frequency is 
mainly determined by the denominator of this equa- 
tion, and below u r by its numerator. For small u 

x"(ku;) w — oj/ikW k ^ 4 Imn(ka;), 

where /ik, Eq. (J2J), is a slowly varying function near 
Q, while cjj^ 4 is strongly peaked at this wave vector 
[see equation (J^}]. Thus, the incommensurate peaks 
may be related only to the last multiplier in this equa- 
tion. The fermion states with the energies £k = ±w/2 
make the main contribution to the peaks in Imll(kaj). 
The location of these states in the momentum space 
is shown in Fig. [3J Here arcs are sections of the hole 
dispersion e-^ at the energies ±0.03t. The blue dashed 
and red dash-dotted lines connect electron-hole pairs, 
which make the main contribution to the magnetic 
susceptibility on the edge and on the diagonal of the 
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Brillouin zone, respectively. The location of these 
states is restricted to the vicinity of the arcs, and the 
difference in the wave vectors of the states is equal to 
the spin-excitation momentum. As a consequence the 
low-frequency susceptibility is sharply peaked at cer- 
tain regions of the Brillouin zone, as seen in Fig. [5J 
In Fig. |3l opposite sides of the magnetic Brillouin 
zone are displaced by the vector Q from each other. 
As seen from the figure, the line segments connecting 
the electron-hole pairs differ from Q and, therefore, 
the response is incommensurate. Notice that there 
are two groups of decay processes for the susceptibil- 
ity on the edge of the Brillouin zone (the two dashed 
lines in Fig. [3]) and only one group for the diagonal 
(the dash-dotted line). This fact explains the larger 
value of x" on the edge in comparison with the diag- 
onal for the parameters of Fig. HJa). As mentioned 
above, in moderately underdoped cuprates the max- 
ima of the low-frequency magnetic susceptibility are 
also located on the zone edges [TU1 EJ H2] ■ 

With increasing oj the influence of the Fermi arcs 
on the susceptibility falls off and x" approaches the 
value it would have in the superconducting phase for 
the same parameters. As a consequence the incom- 
mensurate maxima of the susceptibility have the dis- 
persion of the hourglass shape, which is similar to 
that in the superconducting phase (see, e.g., Fig. 4 
in Ref. [5]). This dispersion is shown in Fig. HJa). 
Notice that this shape of the dispersion is retained 
also for smaller lengths of the Fermi arcs, when the 
susceptibility at very small frequencies is peaked on 
the diagonals [see Fig(2](b)]. As in the superconduct- 
ing phase, the downward-directed branch of the dis- 
persion is connected with the numerator of Eq. (|I1[) . 
and the upward-directed branch with its denomina- 
tor. For uj > ui r the momentum dependence of x" is 
much more isotropic than below the resonance fre- 
quency. Since by definition (fT]) the susceptibility 
coincides with Green's function of spin excitations, 
the upward-directed branch displays the dispersion 
of these excitations. In Fig. [4jb) the frequency de- 
pendencies of %"(kcj) at the maximum and of the 
local susceptibility 

X "(u) = N- 1 J2x"Q<") 
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Figure 4: (a) The dispersion of the maxima in x"(ku;) 
at the edge of the Brillouin zone, k = (k x ,Tr), for 
I = 2 and A = 0.05i. Other parameters are given in 
the text, (b) The value of x"(ko;) at the maximum 
(filled blue circles) and the local susceptibility x"( w ) 
(open red circles) as functions of frequency. 



are shown. As in the superconducting phase, the sus- 
ceptibility maximum is most intensive at the waist 
frequency. The maximum of the local susceptibility 
is shifted to higher frequencies due to a larger sus- 
ceptibility in the central part of the Brillouin zone in 
this frequency range. 

The main difference of the dispersion in Fig. |4] 
from that in the superconducting phase (see Fig. 4 
in Ref. |9 ) is the absence of the superconducting 
gap in the susceptibility at small frequencies. The 
appearance of this gap can be understood in the 
following manner. By analogy with the discussion 
of Fig. [H in the superconducting phase the main 
contribution to the susceptibility at the frequency 
lj < u r is made by fermion states with the energies 
6c = (El + Al) 1 / 2 w oj/2. These states are located 
near the boundaries of crescent pockets, the shape of 
which resembles arcs in Fig. |3l However, in contrast 
to these arcs the length of the pockets shrinks with 
decreasing uj. As a result intensities of the suscepti- 
bility maxima are diminished. Maxima on the edges 
of the Brillouin zone disappear first, since the fermion 
states contributing to them are located on the periph- 
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Figure 5: The susceptibility x"(kiL>) at the incom- 
mensurate maximum on the edge of the Brillouin 
zone as a function of frequency in the pseudogap 
(green circles, I = 2) and superconducting (red 
squares) phases. Blue triangles show the suscep- 
tibility in the pseudogap phase in the case of the 
cusp in the hole dispersion at the Fermi level when 
1 = 2. In all cases A = 0.05i, other parameters are 
given in the text. The experimental susceptibility in 
Lai.875Sro.i25Cu04 measured at T = 4 K [3] is shown 
by purple rhombuses for comparison (the right axis) . 



ery of the pockets. The susceptibility on the diago- 
nals decreases more slowly, as the respective states 
are located near centers of the pockets. Though the 
size of the arcs in Fig. [3] decreases with w, it remains 
finite up to zero frequency. Therefore in contrast to 
the superconducting phase the low-frequency suscep- 
tibility does not vanish in the pseudogap phase. Sus- 
ceptibilities at the incommensurate maxima in the 
two phases as functions of uj are compared in Fig. [5] 

The above-discussed difference in the behavior 
of the low-frequency susceptibility in the supercon- 
ducting and pseudogap phases can account for the 
distinction between the underdoped and overdoped 
cuprates, if additionally a phase separation into su- 
perconducting and pseudogap regions is supposed in 
the former compounds. In contrast to the optimally 
and overdoped crystals, in underdoped cuprates be- 
low T c the superconducting gap in the susceptibil- 
ity is not observed because the low-frequency mag- 
netic response is determined by the pseudogap re- 



gions, while both regions contribute to the suscepti- 
bility for larger uj. Thereby the incommensurability 
parameter - the difference between the wave vector 
of the susceptibility maximum and Q - varies con- 
tinuously from lu = to higher frequencies and the 
dispersion of the maxima has the hourglass shape, as 
in optimally and slightly overdoped cuprates. 

The above discussion allows us to propose the 
mechanism of the reorientation of the low-frequency 
susceptibility maxima from the axial to the diagonal 
direction, which is observed in LSCO for x < 0.05 
|12) . As mentioned previously, for large Fermi arc 
lengths the susceptibility is peaked at the momenta 
(tt, 7r ± k), (it ± k, 7r) up to the smallest considered 
frequencies [see, e.g., Fig. EJa)]. However, the Fermi 
arc length decreases nearly linearly with x [5]. As 
mentioned in the foregoing, the reduction of the arc 
length suppresses the susceptibility on the edge of 
the Brillouin zone more strongly than on the diago- 
nal [see Fig. [2Kb)]. For the used parameters \" on the 
diagonal becomes larger than on the edge for I w 0.6. 
This value approximately corresponds to x = 0.04 [8 , 
which is close to the concentration of the maxima re- 
orientation in LSCO. Notice that in the above calcu- 
lations only the variation in the Fermi arc length was 
taken into account without regard for modifications 
in the hole band with x. The obtained reasonable 
value of the reorientation concentration apparently 
indicates that the Fermi arc length is the main pa- 
rameter which determines this concentration. Let us 
recall that also for such small x the dispersion of the 
susceptibility maxima has the hourglass shape. 

As seen in Fig. [SJ with the used hole dispersion and 
T = the intensity of the incommensurate maximum 
decreases monotonously down to some finite value as 
lo is reduced. In some experiments in underdoped 
cuprates the growth of the susceptibility is observed 
with decreasing oj. As an example of such a behav- 
ior, in Fig. [5] the susceptibility measured in LSCO [2] 
is given. A possible reason for this frequency depen- 
dence may be a small cusp of the hole dispersion on 
the Fermi level. Such a cusp was obtained in the cal- 
culations of the hole dispersion of the 2D t-J model 
for moderate x [13] . A part of the cusp - a bending 
down of the dispersion - is observed in some pho- 
toemission spectra of underdoped cuprates (see, e.g., 
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Figure 6: The section of the hole dispersion along the 
diagonal of the Brillouin zone in the presence (the red 
solid line) and in the absence (the blue dashed line) 
of the cusp. 

Fig. 1 in Ref. [33]). The cusp is shown in Fig. HI It 
increases the combined electron- hole density of states 
at low frequencies in the first term of Eq. ([4]), which 
strengthens the susceptibility. The non-monotonic 
frequency dependence of the susceptibility caused by 
this small change of the hole dispersion is shown by 
filled triangles in Fig. [5] There is a clear similarity 
with the experimental dependence from Ref. [2] . The 
deviation for u> < 1 meV is presumably connected 
with the difference in temperatures in the experiment 
and calculations. 

In Ref. [2] a rapid change of the k width of the 
susceptibility maximum was revealed near the fre- 
quency uj — 4 meV in LSCO with x ~ 0.13 at 
T = 4 K. Results of this work are shown in Fig. [7J by 
green triangles. The widths obtained in our calcula- 
tions for the superconducting and pseudogap phases 
are shown by red circles and blue squares, respec- 
tively. In the superconducting state the susceptibil- 
ity maximum on the edge of the Brillouin zone disap- 
pears below uj — 10 meV. For larger frequencies the 
width grows nearly linearly at first and then, near 
uj r w 0.08t = 40 meV, it starts to decrease, since 
the maxima transform gradually to a sharp resonance 
peak. In the pseudogap phase the number of states 
contributing to the susceptibility maximum is larger 
than in the superconducting state. As a consequence 
the widths of maxima in the former phase are consid- 
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Figure 7: The k width T of the incommensurate 
susceptibility maximum (half width at half maxi- 
mum) as a function of uj in the superconducting 
(red circles) and pseudogap {I = 1.6, blue squares) 
phases. Green triangles show the width measured in 
Lai. 87 Sr . 13 CuO4 at T = 4 K [2 1. The width is given 
in reciprocal lattice units. 



erably larger. As seen in Fig. [JJ there is a steep rise 
of the width near uj = 5 meV. It is connected with 
the shoulder of the main maximum, which moves 
to the corner of the Brillouin zone [see Fig. Ufa)]. 
This shoulder is weak for small uj, however it grows 
rapidly with frequency and determines the width for 
uj > 5 meV. This frequency dependence resembles 
the behavior of the linewidth observed in Ref. [2] . 

As seen in Fig. [U for T = the susceptibility in the 
pseudogap phase tends to some finite value as uj — > 0. 
From Eqs. (j4j and ([TTj) it follows that for finite tem- 
peratures x" starts to decrease for uj\ ~ T and finally 
goes to zero at uj — 0. Usual temperatures of experi- 
ments of quasi-elastic neutron scattering are 2 — 4 K 
(see, e.g., Refs. [2J[3]). Thus, uj 1 ~ 10" 1 meV, which 
is comparable to the energy resolution of these exper- 
iments. Nevertheless, the tendency for a decrease of 
x" as the frequency becomes smaller than T can be 
observed in some experimental spectra, for example, 
in the data reproduced in Fig. [5] 
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4 Conclusion 

Above we supposed that for temperatures T < T c 
the low-frequency magnetic response of undcrdopcd 
cuprates is determined by regions of the pseudogap 
phase in phase-separated crystals. It looks rather rea- 
sonable to extend this picture to the interpretation of 
the influence of modest magnetic fields and impurities 
on magnetic properties of these crystals below T c . An 
external field, which is smaller than the upper critical 
field, enhances the quasi-elastic peaks at incommen- 
surate positions in neutron scattering [3] [22] [23] . In 
optimally doped LSCO such a field leads to a partial 
filling of the superconducting gap in the susceptibil- 
ity [53]. These results suggest that the part of the 
crystal in vertex cores is in the pseudogap phase, and 
in the underdoped case the cores increase the vol- 
ume occupied by this phase. The above supposition 
about the state of the vertex cores was already made 
earlier based on the analysis of tunneling experiments 
[35] [55] and model calculations [37] • Apparently man- 
ifestations of the pseudogap phase are observed not 
only in the vertex cores but also in their vicinity: 
the checkerboard pattern with periodicity of approx- 
imately four lattice spacings was observed inside and 
near the cores by scanning tunneling spectroscopy in 
slightly overdoped Bi 2 Sr 2 CaCu 2 O a; [35] US]- A sim- 
ilar modulation was found above T c , in the pseudo- 
gap phase, in Ref. [30 . As indicated in Ref. [28] . 
the electronic modulation is related to the magnetic 
modulation with twice the wavelength. In this pic- 
ture, the nearly linear magnetic field dependence of 
the quasi-elastic neutron scattering at incommen- 
surate momenta [3] is related to the growing vol- 
ume of the pseudogap regions, which do not overlap, 
while the tendency toward saturation [23] arises with 
their overlapping. The observed strengthening of the 
quasi-elastic scattering by impurities [3T] can be con- 
nected also with regions of the pseudogap phase sur- 
rounding these impurities. 

In summary, the magnetic susceptibility in the 
pseudogap phase of cuprate perovskites was calcu- 
lated using the 2D t-J model and the Mori projec- 
tion operator technique. In contrast to the supercon- 
ducting phase the susceptibility remains finite down 
to very low frequencies, which is the consequence of 



the Fermi arcs in the pseudogap phase. For frequen- 
cies immediately below the resonance frequency the 
magnetic response is incommensurate with maxima 
located on the edge of the Brillouin zone. For even 
lower frequencies the location of maxima depends on 
the length of the Fermi arcs: for long arcs the maxima 
are on the zone edge, while for short arcs they are on 
the diagonals. Together with the experimental fact 
that in LSCO the length of the Fermi arc is propor- 
tional to the hole concentration, this result explains 
the reorientation of the low-frequency susceptibility 
maxima with doping, observed in this crystal. The 
estimated value of the hole concentration, at which 
the reorientation occurs, is close to the experimental 
one. As in the superconducting state, the disper- 
sion of the susceptibility maxima has an hourglass 
shape. Its lower part is connected with the maxima 
in the spin-excitation damping, while the upper part 
reflects the dispersion of spin excitations. The hour- 
glass waist - the resonance frequency - corresponds 
to the gap in the spin-excitation dispersion due to the 
short-range antiferromagnetic order. A small cusp of 
the hole dispersion near the Fermi level, which is con- 
nected with the spin-polaron band, affords a growth 
of the incommensurate susceptibility maxima with 
decreasing frequency. Due to a peculiar shape of the 
susceptibility maximum and its variation with fre- 
quency in the pseudogap phase a sharp change of the 
k width of the maximum occurs at some frequency. 
With the assumption of the phase separation into re- 
gions of superconducting and pseudogap phases in 
underdoped cuprates, these results provide explana- 
tion for the difference in the magnetic response of 
undcrdopcd and overdoped lanthanum and yttrium 
cuprates. Below T c optimally and overdoped crystals 
are characterized by the superconducting gap in the 
magnetic susceptibility, while phase-separated under- 
doped samples have nonzero susceptibility down to 
the lowest measured frequencies at low temperatures. 
Moreover, in the latter crystals the susceptibility may 
grow with decreasing frequency. These peculiarities, 
as well as the observed rapid change of the k width 
of the susceptibility maximum with frequency are re- 
produced in reasonable agreement with experiment, 
if we relate them to the pseudogap regions of the 
phase-separated crystal. 



9 



This work was supported by the project 
TLOFY0145 and by the European Regional 
Development Fund (project TK114). 

References 

[I] J. Rossat-Mignot, L.P. Regnault, C. Vettier, 
P. Bourges, P. Burlet, J. Bossy, J.Y. Henry, 

G. Lapertot, Physica C 185-189, 86 (1991) 

[2] M. Kofu, S.-H. Lee, M. Fujita, H.-J. Kang, 

H. Eisaki, K. Yamada, Phys. Rev. Lett. 102, 
047001 (2009) 

[3] D. Haug, V. Hinkov, A. Suchaneck, D.S. Inosov, 
N.B. Christensen, Ch. Niedermayer, P. Bourges, 
Y. Sidis, J.T. Park, A. Ivanov, C.T. Lin, J. Mesot, 
B. Keimer, Phys. Rev. Lett. 103, 017001 (2009) 

[4] M. Lavagna, G. Stemmann, Phys. Rev. B 49, 
4235 (1994) 

[5] J. Brinckmann, P.A. Lee, Phys. Rev. Lett. 82, 
2915 (1999) 

[6] A. Sherman, M. Schreiber, Phys. Rev. B 68, 
094519 (2003) 

[7] M.R. Norman, H. Ding, M. Randeria, J.C. Cam- 
puzano, T. Yokoya, T. Takeuchi, T. Taka- 
hashi, T. Mochiku, K. Kadowaki, P. Guptasarma, 
D.G. Hinks, Nature 392, 157 (1998) 

[8] T. Yoshida, X.J. Zhou, D.H. Lu, S. Komiya, 
Y. Ando, H. Eisaki, T Kakeshita, S. Uchida, 
Z. Hussain, Z.-X. Shen, A. Fujimori, J. Phys.: Con- 
dens. Matter 19, 125209 (2007) 

[9] A. Sherman. larXiv:1108.4179l (2011) 

[10] H. Yoshizawa, S. Mitsuda, H. Kitazawa, K. Kat- 
sumata, J. Phys. Soc. Jpn. 57, 3686 (1988) 

[II] R.J. Birgeneau, Y. Endoh, Y. Hidaka, K. Kaku- 
rai, M.A. Kastner, T. Murakami, G. Shirane, 
T.R. Thurston, K. Yamada, Phys. Rev. B 39, 2868 
(1989) 



[12] M. Fujita, H. Hiraka, M. Matsuda, M. Matsuura, 
J.M. Tranquada, S. Wakimoto, G. Xu, K. Yamada, 
larXiv:1108.443Tl (2011) 

[13] A. Sherman, M. Schreiber, Phys. Rev. B 55, 
R712 (1997) 

[14] K. Terashima, H. Matsui, T. Sato, T. Takahashi, 
M. Kofu, K. Hirota, Phys. Rev. Lett. 99, 017003 
(2007) 

[15] H. Mori, Progr. Theor. Phys. 34, 399 (1965) 

[16] J. Kondo, K. Yamaji, Progr. Theor. Phys. 47, 
807 (1972) 

[17] H. Shimahara, S. Takada, J. Phys. Soc. Japan 
61, 989 (1992) 

[18] A. Sherman, M. Schreiber, Eur. Phys. J. B 32, 
203 (2003) 

[19] A.K. McMahan, J.F. Annctt, R.M. Martin, 
Phys. Rev. B 42, 6268 (1990) 

[20] V.A. Gavrichkov, S.G. Ovchinnikov, 
AA. Borisov, E.G. Goryachev, JETP (Rus- 
sia) 91, 369 (2000) 

[21] A. Sherman, M. Schreiber, in: Pseudogap in- 
High Temperature Superconductors - Studies of 
High Temperature Superconductors, vol. 27, edited 
by A.V. Narlikar (Nova Science Publishers, New 
York, 1999), p. 163 

[22] S. Katano, M. Sato, K. Yamada, T. Suzuki, 
T. Fukase, Phys. Rev. B 62, R14677 (2000) 

[23] B. Lake, H.M. R0nnow, N.B. Christensen, 
G. Aeppli, K. Lefmann, D.F. McMorrow, 
P. Vorderwisch, P. Smeibidl, N. Mangkorntong, 
T. Sasagawa, M. Nohara, H. Takagi, T.E. Mason, 
Nature 415, 299 (2002) 

[24] B. Lake, G. Aeppli, K.N. Clausen, D.F. McMor- 
row, K. Lefmann, N.E. Hussey, N. Mangkorntong, 
M. Nohara, H. Takagi, T.E. Mason, A. Schroder, 
Science 291, 1759 (2001) 



10 



[25] Ch. Renner, B. Revaz, K. Kadowaki, I. Maggio- 
Aprile, 0. Fischer, Phys. Rev. Lett. 80, 3606 
(1998) 

[26] G. Levy, M. Kugler, A. A. Manuel, 0. Fischer, 
Phys. Rev. Lett. 95, 257005 (2005) 

[27] C. Berthod, B. Giovannini, Phys. Rev. Lett. 87, 
277002 (2001) 

[28] L.E. Hoffman, E.W. Hudson, K.M. Lang, 
V. Madhavan, H. Eisaki, S. Uchida, J.C. Davis, 
Science 295, 466 (2002) 

[29] K. Matsuba, S. Yoshizawa, Y. Mochizuki, 
T. Mochiku, K. Hirata, N. Nishida, J. Phys. Soc. 
Jpn. 76, 063704 (2007) 

[30] M. Vershinin, S. Misra, S. Ono, Y. Abe, 
Y. Ando, A. Yazdani, Science 303, 1995 (2004) 

[31] M. Fujita, M. Enoki, S. Iikubo, K. Kudo, 
N. Kobayashi, K. Yamada. la7Xiv:0903.539"i1 ([2009) 



11 



